Equilibrium Circular Photogalvanic Effect in a hybrid superconductor-semiconductor 

system 
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A dc electric current can be induced in a hybrid semiconductor-superconductor system under 
illumination it by a circularly polarized light with the frequency below the energy of semiconductor 
interband transitions. In conditions when the light beam is unable to create real electron-hole 
excitations, this phenomenon is reminiscent of the Meissner effect in the static magnetic field. Such 
an effect can be employed in systems combining cavity photons and superconducting quantum 
circuits. 

PACS numbers: 72.40.+W, 74.45.+C, 75.70.Tj 



In noncentrosymmetric semiconductors a circularly po- 
larized light can generate the dc electric current This, 
so called, circular photogalvanic effect (CPE) is deter- 
mined by the second-order response to the electric field 
of the incident light. It involves spin polarization of elec- 
trons excited from valence bands split by the spin-orbit 
interaction (SOI) into bands with the total angular mo- 
ments 3/2 and 1/2. The polarized carriers, in their turn, 
give rise to the electric current, due to SOI associated 
with the lack of inversion symmetry. This phenomenon 
attracted much attention recently [2J in connection with 
perspective spintronic applications. 

As expected, in the case of a dissipative transport in 
normal electron systems at stationary conditions, CPE 
takes place only if light illumination creates a thermo- 
dynamically nonequilibrium distribution of carriers in a 
semiconductor On this reason, one can not expect 
this effect to occur in bulk semiconductors, if the fre- 
quency of the incident light is less than the energy gap 
between the valence and conduction bands and when in- 
direct phonon/impurity assisted optical transitions are 
weak. On the other hand, in superconductors the flow 
of electrons might be created without exciting the sys- 
tem from its thermal equilibrium. Such a phenomenon 
becomes possible due to a macroscopic coherency of the 
condensate wave function, so that spacial variations of 
its phase determine the supercurrent. The well known 
example is the Meissner effect, where the static magnetic 
field can not produce electron-hole excitations. It, how- 
ever, induces a superconducting electric current. The 
goal of this Letter is to show that CPE can be observed 
in a hybrid superconductor-semiconductor system even if 
the incident light does not drive it from the thermal equi- 
librium. The electric current is induced due to changes 
in the spectrum and wave function of the many-electron 
system, rather than from deviation of the electron's dis- 
tribution function from its thermal equilibrium. Apart 
from the fundamental interest, this phenomenon has a 
practical value when it is employed in nanoscale opto- 
electronic devices, because it allows to reduce dissipation 
losses compare to CPE in normal systems. It can also be 



an efficient tool to combine microwave electrodynamics in 
superconducting quantum circuits [i [with cavity optical 
fields coupled to atomic ensembles 5j, like the recently 
observed [6| interaction of spin ensembles to supercon- 
ductor resonators. This suggests evident connections to 
many topics of current interest, ranging from atomic and 
polaritonic Bose condensates to spintronics and quantum 
information processing. 

The equilibrium CPE takes place in a hybrid 
semiconductor-superconductor system due to the prox- 
imity effect which induces Cooper pair correlations in 
the semiconductor. Various hybrid systems displaying 
the strong proximity effect has been recently fabricated 
[7|. In order to demonstrate main features of CPE and 
make our quantitative analysis more transparent, a sim- 
ple model will be considered here, where a noncentrosym- 
metric semiconductor film is in a planar contact with su- 
perconductor having the singlet order parameter. The 
film contains an n-doped quantum well (QW) close to 
the metal surface (Fig. 1). Note, that such a sandwich 
system is now being considered as a key element of a 
topological superconductor [8(. 

Let us assume that the circularly polarized light with 




FIG. 1: (Colour online) A sketch of the system. An optically 
transparent semiconductor film (top) is in a contact with a 
superconductor (bottom). A thin layer between them depicts 
a doped quantum well. An incident electromagnetic wave 
induces in this well an electric current J perpendicular to the 
vector P~ExE', where E is the electromagnetic field. The 
light beam can be also incident from the metal side, if the 
metal film is thinner than the skin-layer. 
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the frequency Ui is incident onto the semiconductor sur- 
face. The vector potential of the total electromagnetic 
field in the QW region is represented by its Fourier 
components A and A*, corresponding to frequencies 
u) = u>i and u> — —0Ji, respectively. It will be assumed 
that the electromagnetic field causes only direct tran- 
sitions between valence and conduction bands. Hence, 
phonon and impurity assisted indirect transitions, as 
well as a very small momentum transfer by photons, 
will be ignored. The corresponding matrix elements 
are given by M av {uji) = (e/mc)(k, a|p • A|k, v) and 
M av (—u!i) = (e/mc)(k, a\p ■ A*|k,u) , where k denotes 
the two-dimensional wave vector of QW electrons and p 
is their momentum operator, a is the conduction elec- 
tron spin variable and v labels the valence bands. In 3D 
superconductors the latter are the heavy-hole and light- 
hole bands with the total angular moment I = 3/2 and its 
projections ±3/2 and ±1/2, respectively, plus the split- 
off band with I = 1/2, l z = ±1/2. In QW the label v 
runs also through valence subband indexes, while only 
the lowest conduction subband is assumed to be occu- 
pied. In noncentrosymmetric semiconductors SOI splits 
the energies of bands with opposite angular moments. 
For CPE this splitting is crucial. In the case of a n- 
doped QW the conduction band splitting appears to be 
most important, although at strong resonance conditions 
the hole splitting can be equally important. Therefore, 
the latter will be ignored under the assumption that it 
is much less than the resonance detuning. At the same 
time, the conduction band splitting is determined by the 
SOI Hamiltonian H so = <r ■ hk, where a = (a x ,a y ,a z ) 
is the vector of Pauli matrices and hk = — h k is the 
spin-orbit field. The latter will be assumed to have only 
x and y components. This takes place for the Rashba 
field, as well as for the Dresselhaus field in [001] oriented 
QW |. 

The semiconductor film contacts 3D metal through a 
high enough tunneling barrier, so that the broadening of 
QW states due to resonance with 3D continuum is much 
smaller than the metal superconducting gap. Therefore, 
in the leading approximation, the electric current in QW 
can be calculated by ignoring a leaking of electrons into 
the metal. The system is assumed to be clean enough, 
with the mean electron scattering rate in QW to be much 
smaller than the proximity induced energy gap and h^. 
Hence, the Green's functions are diagonal with respect 
to electron wave vectors. Considering a thermodynamic 
equilibrium with the temperature T, the stationary elec- 
tric current density can be expressed in terms of the equi- 
librium Keldysh function, as 
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where G^\oj) are the retarded and advanced Green's 
functions. The index n labels the energy bands including 
conduction n ~ c and valence n = v bands. Although 
the valence bands are far from the Fermi level, their in- 
clusion into the current expression is necessary, because 
they contribute to important spectral corrections asso- 
ciated with the coupling of electrons to the electromag- 
netic field. The current operator of conduction electrons 
is given by 
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while the valence band currents are given by correspond- 
ing band velocities de v .^/dk, where e Mi k are the band 
energies measured with respect to the chemical poten- 
tial. The above expressions are written in the Gor'kov- 
Nambu representation, where the new electron destruc- 
tion operators Ck^^i = c^ n>(7 and c k: „,<T : 2 = C-u,n,s are 
introduced, with Pauli matrices ti,T2,T3 acting in the 
space of indices 1 and 2. As mentioned in Ref.[10|. this 
basis is more convenient for studying the spin-dependent 
transport in superconducting systems. 

The next step is to specify important contributions 
to the electron self-energy due to the superconducting 
proximity effect and interaction with the electromagnetic 
field. The proximity effect is represented by a nondiag- 
onal in the Nambu space contribution to the self-energy 
of QW electrons. In the case when parallel to the in- 
terface wave vectors of tunneling particles are conserved, 
this self-energy has the form 
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where Xk,fc z is the tunneling matrix element, with k z de- 
noting the normal component of the wave vector in the 
metal side of the interface. Assuming the real order pa- 
rameter Aq, the Green's function in ((3]) can be written 



as G 
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(ui — T^ek — a z TiAo ± iS) . For a wide-band 
metal, where Tk.fc, slowly varies near the Fermi level, the 
self-energy is expressed as 
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where T^(e) = J2 k iTk^J 2 ^ — e^.) determines a finite 
confinement lifetime of QW electrons due to resonance 
with the continuum of metal states. The proximity ef- 
fect is represented by the second term in the numerator 
of Eq. @ . Since we assumed that Tk <C Ao and the tem- 
perature is low, one can neglect ui in comparison with Ao. 
In this case the gap A in the spectrum of QW electrons is 
given by A = Tk(0) and it is reasonable to neglect a very 
weak dependence on k, because the Fermi wave vector is 
much smaller in QW than in the metal. 

The stationary effect of the electromagnetic field onto 
the electron current is determined by the second-order 
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contribution to the conduction electron self-energy. Its 
11 component in the Nambu space is 

£ii, oj8 (w J k) = 53[M au (w i )M| (wi)G ! Sk(w + w i ) + 

McoC-WijM^f-WjJGSktw - ut)] , (5) 

where the unperturbed valence-band Green's functions 
are G^ k (u;) = (v — e^k) -1 - From the above defini- 
tion of Gor'kov-Nambu operators, £22 can be expressed 
as £22,0/3 (k- 1 , k) = — £ xl s & {— cj, — k). The valence-band 
wave functions who determine the matrix elements in Eq. 
(|5|) are linear combinations of the functions | J)|ct), where 
J = X,Y, or Z denote / = 1 orbitals and a is the spin 
index. Therefore, depending on the spin variables a and 
/3 in Eq. ((SJ), the electromagnetic field enters as various 
combinations of A x A*, A • A* and A Z A* Z . The circular 
photogalvanic effect is associated with P = iA x A*. 
In general, its contribution to Eq. has the form 

1] CP • er, where the factor C contains resonance de- 
nominators. The resonance detuning, however, was as- 
sumed to be much larger than A and h^. Since U) ~ 
max(A, /ik, IibT), a weak dependence of the self-energy 
on u> can be ignored and £ becomes 



£(w,k) = r 3 (P-<x)C k . 



(6) 



The electric current density will be calculated in the 
leading order with respect to £(w,k). Hence, the corre- 
sponding correction to the Green's function of conduction 
electrons in Eq. (JTJ can be written as 

5G r k {a \u) = G£ (o) (w)E(w,k)G£ (o) M , (7) 

The unperturbed functions, in their turn, are given by 

G° r(a) H = (u>-ne ck -cr-h k -Y,;( a \u;,k)±i5)- 1 , (8) 

where, as was mentioned above, £^ a '(a>,k) ~ — ri<r z A. 

The Green's functions of valence electrons in Eq. (fTJ) 
are calculated in a similar way. The electromagnetic field 
contributes to their self-energy through the same matrix 
elements as in Eq. ([3]), with the intermediate states com- 
ing from the conduction band. It is easy to show that 



Tr [j„ k G„kM] = Tr 
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Let us choose the x- axis parallel to the xy projection 
of P. For simplicity, the conduction and valence bands 
are assumed to be isotropic and SOI is taken in the form 
of the Rashba interaction, where h x = jky^hy — —"fk x . 
Substituting Eq.© and Eq. into Eq. after some 
algebra one arrives to J x = and 
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where E = Ve 2 + A 2 , 6 = (1 + dh\C k / dh\k) k=kF and 
h = "fkp. As expected, the equilibrium CPE turns to 
in the normal state |3j]. One can immediately see it 
from Eq. pU|) at A — > 0. In the opposite limit of the 
large superconducting gap A 3> k B T the current reaches 
its maximum magnitude. In this case cosh~ 2 (E/k B T) in 
the integrand is exponentially small and in the leading 
approximation J y is given by 



J y — 



eN Fl C kF P x f( 2 _^l 
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where £ = y/l + (A 2 /h 2 ). We note that the expression 
in brackets smoothly varies between 2b— 2, at h <C A and 
26, at h ^> A. Therefore, if the light frequency is not too 
close to the inter-band transition energy, one gets 6 ~ 1 
and this expression is of the order of unity. 

It is important to note that CPE is determined by the 
dependence of the electron-photon second-order scatter- 
ing on the electron spins. The corresponding self-energy 
£(w, k) results in the spin orientation of electrons in the 
direction parallel to P and, being combined with SOI, 
gives rise to the electric current. A basic distinction of 
the equilibrium CPE in superconductors is that at low 
temperatures there are no single-particle spins to be po- 
larized. Their role, however, is played by triplet Cooper 
pairs, which admix to singlets due to SOI. This situation 
resembles the effect of the Zeeman interaction, which to- 
gether with SOI also leads to the electric current in su- 
perconducting systems fill ]. 

Some comments are needed regarding the vector P in 
Eqs. (|6I10|) . In the considered geometry (Fig.l), the 
electromagnetic field near the QW is a sum of the in- 
cident and reflected waves. Therefore, in the presence 
of highly reflecting metal, this vector is strongly modi- 
fied with respect to its plane-wave value. Let the inci- 
dent and reflected light beams in the x'z plane. Near 
QW the circularly polarized incident wave has the elec- 
tric field components E\ = — sin OEq , E x , — cos 0E$ and 
EL = iEq, where 9 is the incident angle at the metal- 
semiconductor interface. In general, the total electric 
field E near the interface may be expressed in terms of 
the surface imp edance Z(uj.i) . From these well known 
expressions [12j, the in-plane components of the vector 
P = i(c 2 /cj 2 )(E x E*) may be easy calculated: 



P x > 



LEol 2 — .rsinflcos* 9 Re 
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where the terms Z cos 9 have been neglected, taking into 
account that the impedance of a highly reflective and 
thick enough metal film is small. In the range of not too 
small frequencies, Z is represented by its Fresnel expres- 
sion Z(u>i) = y/ ' —iuiiCs (u>i )/4ira(uii), where cr(wj) and 
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e s (uJi) are the metal conductivity and dielectric function 
of the semiconductor film, respectively. At oji ^> A the 
former is almost the same as that of the normal metal. 
Both, £ s (uii) and cr(wi) depend on the frequency and 
many other factors. Therefore, it is useful to consider 
typical situations. In the case of a highly reflective metal 
and e s {cui) > we have Re(Z) < Im(Z). Hence, Eq. (fT2|) 
gives the vector P perpendicular to the x'z plane. Since 
the x-axis in Eq. (fT0|) has been chosen parallel to P, this 
means that the electric current will be directed parallel 
to x'-axis. At the same time, for relatively low reflective 
metals the real and imaginary parts of Z will be of the 
same order of magnitude. Therefore, the vector P and, 
hence, the electric current will be directed arbitrary in 
the x'y' plane, their directions varying with the light fre- 
quency. It should be noted that a setup different from 
Fig. 1 might be chosen. For example, the light beam 
can be incident from the superconductor side, penetrat- 
ing through a semi-transparent film whose thickness is 
less than the skin-layer depth. In this case the vector P 
will have quite different from Eq. ([T2j) components. 

For an order-of-magnitude evaluation of CPE, the fac- 
tor C may be approximated by its bulk expression [lj, 
that is valid far enough from the resonance: 

where E g is the semiconductor energy gap and A so de- 
notes the split-off energy. Taking u)\ ~ E g ~ A so ~ 
0.4 eV (InAs), we get C ~ (e 2 /mc 2 )(p 2 cv /mE g ) ~ 
(e /mc 2 )(m/m*), where m* is the effective electron 
mass. The Rashba parameter 7 can vary depending on 
QW characteristics. For InAs it can be more than 10~ n 
eVm We take 7 =10~ 12 cVm, fixing thus h around 
lmeV at kp = 10 6 cm -1 , that is much less than the Fermi 
energy of QW electrons. At P x ~ c 2 \Ei\ 2 /uj 2 and the 
moderate electric field strength Ei = 10 4 V/cm, Eq.(fTTj) 
gives J y ~ 1/iA/cm. This value will increase considerably 
at resonance conditions, by the factor E g /5 res , where 
5 res is a detuning from the resonance. The above theory 
restricts 5 res by a much larger value than the electron 
and hole spin-orbit splittings. In the safe range of rela- 
tively large detunings ~10meV, the resonance enhance- 
ment factor is about 40 for InAs based QW. 

A really dramatic enhancement can be reached, if the 
incident light intensity is periodically modulated with a 
microwave frequency that is close to an eigenfrequency 
of a quantum circuit incorporating the considered super- 
conducting system. In principle, such a modulation could 
be provided by the Rabi splitting of an optical cavity 
mode strongly coupled to atomic ensembles having spin 
resolved electron transitions. In this way CPE gives rise 
to a nonlinear interaction of optical cavity fields + asso- 



ciated atomic ensembles with quantum circuits. 

In conclusion, the electric current induced by the cir- 
cular photogalvanic effect has been calculated in the case 
of a superconducting system. It has been shown that, 
unlike CPE in a normal system, this effect becomes pos- 
sible without driving electrons out of the thermodynamic 
equilibrium and without mediation of the impurity or 
phonon scattering, as well as intersubband optical tran- 
sitions. This phenomenon has been considered for a 
noncentrosymmetric semiconductor QW, where super- 
conductivity is induced by the proximity effect and the 
spin-orbit interaction is represented by the Rashba term. 
In the considered setup the magnitude and direction of 
the electric current can be varied by changing the inci- 
dent light frequency and polarization. 
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